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Abstract. We describe a method for computing the highest de- 
gree coefficients of a weighted Ehrhart quasi-polynomial for a ra- 
tional simple polytope. 



1. Introduction 

Let p be a rational polytope in V = M*^ and h{x) a polynomial 
function on V. A classical problem is to compute the sum of values of 
h{x) over the set of integral points of p, 

S{p,h)= J2 h{x). 

The function h{x) is called the weight. When p is dilated by an integer 
n G N, we obtain a function of n which is quasi-polynomial, the so- 
called weighted Ehrhart quasi-polynomial of the pair (p, h) 



S{np,h) = J2E„ 



d+M 

It has degree d+M, where = deg h. The coefficients are periodic 
functions of n G N, with period the smallest integer q such that qp is a 
lattice polytope. 

In [1], Barvinok obtained a formula relating the kth coefficient of 
the (unweighted) Ehrhart quasi-polynomial of a rational polytope to 
volumes of sections of the polytope by affine lattice subspaces parallel 
to k-dimensional faces of the polytope. As a consequence, he proved 
that the fco highest degree coefficients of the unweighted Ehrhart quasi- 
polynomial of a rational simplex can be computed by a polynomial 
algorithm, when the dimension d is part of the input, but kg is fixed. 
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The sum S{p,h) has natural generahzations, the intermediate sums 
S^{p,h), where L C ^ is a rational vector subspace. For a polytope 
p G V and a polynomial h{x) 

s'^{p,h) = J2 [ h{y)dy, 

^ Jpn{x+L) 

where the summation index x runs over the projected lattice in V/L. 
In other words, the polytope p is sliced along lattice affine subspaces 
parallel to L and the integrals of h over the slices are added up. For 
L = V, there is only one term and S^{p, h) is just the integral of h{x) 
over p, while, for L = {0}, we recover S{p, h). Barvinok's method was 
to introduce particular linear combinations of the intermediate sums, 

Y,p{L)S\p,h). 

Lec 

It is natural to replace the polynomial weight h{x) with an exponential 
function x ^— e^^'^\ and consider the corresponding holomorphic func- 
tions of ^ in the dual V* . Moreover, one can allow p to be unbounded, 
then the sums 

5^(p)(0 = E / 

J, Jpr\{x+L) 

still make sense as meromorphic functions on V* . The map p ^ 
S^iV^iC) is a valuation. In [2j, we proved that Barvinok's valuation 
Slg£ P(-^)'S'^(p)(0 approximates S'(p)(^) in a sense which is made 
precise below. As a consequence, we recovered Barvinok's main theo- 
rem of [4j and we sketched a method for computing the highest degree 
coefficients of the Ehrhart quasipolynomial of a rational simplex which 
is hopefully easier to implement. The proof in [2] relied on our Euler- 
Maclaurin expansion of these functions. 

The main interest of the present article is to give a simpler formula- 
tion and an elementary proof of the approximation result of [2] , in the 
case of a simple polytope. 

In a forthcoming article, in collaboration with J. De Loera and M. 
Koeppe, we plan to apply the results of this article to derive a poly- 
nomial algorithm for computing the highest degree coefficients of a 
weighted Ehrhart quasi-polynomial relative to a simplex, when and 
the degree of the weight h{x) are fixed, but the dimension of the simplex 
is part of the input. The article [1] dealt with the integral h{x)dx, 
which is of course the highest coefficient, if h{x) is homogeneous. 

To explain the main idea, let us assume in this introduction that 
the vertices s of p are integral. Using a theorem of Brion, one writes 
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S'(p)(^) as a sum of the generating functions of the supporting cones 
at the vertices s of p, 

sev(p) 

Then the dilated polytope np has vertices ns with the same cone of 
feasible directions c^, thus 

sev{p) 

The generating function of a cone c is a meromorphic function of a 
particular type, namely, near ^ = 0, it is a quotient of a holomorphic 
function by a product of linear forms. Hence, it admits a decomposition 
into homogeneous components. One shows that the lowest degree is 
> -d. 

Let us fix an integer , < ko < d. Our goal is to compute the k^+l 
highest degree coefficients of the Ehrhart quasi-polynomial of p for the 
weight h{x) = . As we explain in Section 3, this computation 

amounts to computing the lowest homogeneous components 

5'(c,)[_d+fc](0, k = 0,...,ko, 

of the generating functions of the cones c^. 

Our main result. Theorem [121 is an expression, depending on ko, 
for the components S{s + c)[^d+k]{0^ = 0, . . . , fco in the particular 
case where the cone s + c is simplicial. For a unimodular cone, the 
generating function is given by a "short" formula, thus its lowest degree 
components are readily computed. In general, let Vi G Z'^, i = 1, . . . ,d, 
be integral generators of the edges of c. The finite sum 

/(o = E ^^^'^^ 

is an analytic function of ^. (If c is unimodular and the generators Vi 
are primitive, then /(^) = 1). The term of degree k of /(^) is given by 

^ k\ ' 

^6(Eti[o,i[f»)nzd 

A monomial of degree k in can involve at most k vari- 

ables among the ^j. From this elementary remark, we deduce that 
the terms of degree < ko of /(^) can be computed using a Moebius 
type combination of sums similar to f{C), in dimension < /cq, and de- 
terminants. As a consequence, we obtain an expression for the terms 



4 VELLEDA BALDONI, NICOLE BERLINE, AND MICHELE VERGNE 

S{s + c)[_rf_|_fc](^), k — 0, . . . ,ko, which involves the generating func- 
tions of cones in dimension < kg only. This feature is useful because, 
when the dimension is fixed, there is an efficient algorithm for decom- 
posing a simplicial cone into a signed combination of unimodular cones, 
due to Barvinok. 



2. Notations and basic facts 

2.1. We consider a rational vector space V of dimension d, that is to 

say a finite dimensional real vector space with a lattice denoted by A. 
We will need to consider subspaces and quotient spaces of V, this is 
why we cannot just let = R"^ and A = Z'^. The set A ® Q of rational 
points in V is denoted by Vq. A subspace L of F is called rational if 
L n A is a lattice in L. If L is a rational subspace, the image of A in 
V/L is a lattice in V/L, so that V/L is a rational vector space. The 
image of A in V/L is called the projected lattice. 

A rational space V, with lattice A, has a canonical Lebesgue measure 
dx, for which V/A has measure 1. 

A convex rational polyhedron p in (we will simply say polyhe- 
dron) is, by definition, the intersection of a finite number of half spaces 
bounded by rational affine hypcrplancs. We say that p is solid (in V) 
if the affine span of p is 

In this article, a cone is a polyhedral cone (with vertex 0) and an 
affine cone is a translated set s + c of a cone c. 

A polytope p is a compact polyhedron. The set of vertices of p is 
denoted by V(p). For each vertex s, the cone of feasible directions at s 
is denoted by c^. 

A cone c is called simplicial if it is generated by independent elements 
of V. A simplicial cone c is called unimodular if it is generated by 
independent integral vectors vi, . . . ,Vk such that {vi, . . . ,Vk} can be 
completed to an integral basis of V. An affine cone o is called simplicial 
(resp. simplicial unimodular) if it is the translate of a simplicial (resp. 
simplicial unimodular) cone. 

2.2. Generating functions. 

Definition 1. We denote by T-C{V*) the ring of holomorphic functions 
defined around E V* . We denote by Ai{V*) the ring of meromorphic 
functions defined around E V* and by M.i{V*) C M.{V*) the subring 
consisting of those meromorphic functions (f){^) such that there exists 
a product of linear forms D{^) such that D(^)0(^) is holomorphic. 
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A function 0(^) G AieiV*) has a unique expansion into homogeneous 
rational functions 

m>'mo 

If P is a homogeneous polynomial on V* of degree p, and D a product 
of r linear forms, then ^ is an element in Aii{V*) homogeneous of 
degree m = p — r. For instance, |^ is homogeneous of degree 0. On this 
example, we observe that a function in Ai^iy*) which has no negative 
degree terms need not be analytic. 

Let us recall the definition of the functions /(p) and S{p) G Aii{V*) 
associated to a polyhedron p, (see for instance the survey [5]). 

Proposition 2. There exists a unique map I which to every polyhedron 
p G V associates a meromorphic function /(p) G M-iiV*), so that the 
following properties hold: 

(a) If p is not solid or if p contains a straight line, then l{p)=0. 

(b) If ^ E V* is such that e^^'^^ is integrable overp, then 

lipm = [ e^^^^^dx. 
Jp 

(c) For every point s G Vq, one has 

J(. + p)(0 = e«'^>/(p)(0. 

Proposition 3. There exists a unique map S which to every polyhedron 
p (ZV associates a meromorphic function S{p) G M.^{y*), so that the 
following properties hold: 

(a) If p contains a straight line, then S{p)=0. 

(b) If ^ ^ V* is such that e^^'^^ is summable over the set of lattice 
points ofp, then 

s{pm = E 

xe pnA 

(c) For every point s G A, one has 

sis+pm = e^^''^s{pm- 

Moreover the maps p t-^ /(p) and p i— > ^(p) have additivity proper- 
ties, with consequence the fundamental Brion's theorem. 

Theorem 4 (Brion, [7]). Let p be a polyhedron with set of vertices 
V(p). For each vertex s, let Cs be the cone of feasible directions at s. 
Then 

S{p) = E S{s + Cs) andl{p)= ^ ^(-s + c,). 
sev(p) sevip) 
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2.3. Notations and basic facts in the case of a simplicial cone. 

Let Vi,i — l,...,d be linearly independent integral vectors and let 
c = Yli=i '^'^Vi be the cone they span. 

Definition 5. The semi-closed unit cell B of the cone (with respect to 
the generators Vi,i — 1, . . . ,d) is the set 

d 

B = ^[0,lbi. 

i=l 

We recall the following elementary but crucial lemma. 

Lemma 6. (i) The affine cone (s + c) fl A is the disjoint union of the 

translated cells s + B + f , for v G X]j=i ^"^i- 

(a) The set of lattice points in the affine cone s + c is the disjoint union 
of the sets x + Yl'i=i when x runs over the set (s + B) fl A. 
(Hi) The number of lattice points in the cell s + B is equal to the volume 
of the cell with respect to the Lebesgue measure defined by the lattice, 
that is 

Card((s + B) n A) = \det{vi)\. 



Let s eVq. We have immediately 



1Y\ det{vi 



(1) Hs + cm-e 

The study of the function S{s + c)({) will be the main point of this 
article. It reduces to the study of the holomorphic function 5'(s-|-B)(^) 
defined by the following finite sum, over the lattice points of the unit 
cell. 

Definition 7. 

a;e(s+B) nA 

Lemma 8. 

(2) S{s + cm = Sis + Bm , / . 

11^=1(1 - 

In particular, S{s + c) G M.i{V*), thus it admits a decomposition into 

homogeneous components, 

(3) S{s + c) (0 = S[_^ (s + c) (0 + ^[-d+i] {s + cm + ..., 
and the lowest degree term S[^j^{s + c)(^) is equal to /(c) (^) 
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Proof. (I2l) follows from Lemma [6] (ii). Next, we write 



(4) Ut^ = IIt 



The function is holomorphic with value —1 for a; = 0. Thus 
S{s + c) G J^i(y*). The value at ^ = of the sum over the cell is the 
number of lattice points of the cell, that is the volume | det(fj)|. This 
proves the last assertion. □ 

3. Weighted Ehrhart quasipolynomials 

Let p C y be a rational polytope and let h{x) be a polynomial 
function of degree M on V. We consider the following weighted sum 
over the set of lattice points of p, 

E ^(^)- 

xepnA 

When p is dilated by a non negative integer n, we obtain the quasi- 
polynomial of the pair (p, h), 

d+M 

E h{x)= Y.EmU"'. 
xgnpnA m=0 

The coefficients Em actually depend on n, but they depend only on n 
mod q, where q is the smallest integer such that gp is a lattice polytope. 
If h{x) is homogeneous of degree M, the highest degree coefficient E^+m 
is equal to the integral h{x)dx. 

Let us fix a number ko. Our goal is to compute the fco + 1 highest 
degree coefficients E^, for m = M + d, . . . , M + d — k^. 

We concentrate on the special case where the polynomial h{x) is a 
power of a linear form 

n{x) = ; — . 

Of course, any polynomial can be written as a linear combination of 
powers of linear forms. 

We will explain our results with the simplifying assumption that the 
vertices of the polytope are lattice points. 

Definition 9. We define the coefficients -Em(p, ^, M), m = 0, . . . , M+d 

by 



2^ Ml ^ 

xGnpnA m=0 



E^(p,e,M)n™. 
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Proposition 10. Let p he a lattice polytope. Then, for k > 0, we have 

{M + d-k)\' 



If „\M+d-k 

(5) i?M+,_,(p,e,M)= J] 



The highest degree coefficient is just the integral 



Remark 11. As functions ofC,, the coefficients Em{p,C,,M) are poly- 
nomial, homogeneous of degree M. However, in they are expressed 
as linear combinations of rational functions of C,, whose poles cancel 
out. 

Proof. The starting point is Brion's formula. As the vertices are lattice 
points, we have 

(6) E ^^^'^^= E S{s + cM)= E ^^^''^S{cM)- 

xepnA seV(p) sgV(p) 

When p is replaced with np, the vertex s is replaced with ns but the 
cone Cs does not change. We obtain 

E ^^'"^ = E ^""^'''^sicM)- 

We replace ^ with tC,, 

x&ipnA seV(p) 
The decomposition into homogeneous components gives 

sicsM) = rt(c.)(0 + t"'+'5M+i](c.)(0 + ■ ■ ■ + t'S[k]{csm + ■■■■ 

Hence, the t^-term in the right-hand side is equal to 
Thus we have 

^' ^ Ml ^ {M + d)\^''^^' 

xenpnA seV(p) ^ ' 



If: „\Af+ci-l 
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On this relation, we read immediately that ^^^.g^pnA mi ^ polyno- 
mial function of n of degree M + and that the coefficient of n^+'^-^ 
is given by ([5]). The highest degree coefficient is given by 

sev(p) ^ ^' 

Applying Brion's formula for the integral, this is equal to the term of 
^-degree M in /(p)(^), which is indeed the integral dx. □ 

From Proposition [TDl we draw an important consequence: in order 
to compute the ko + 1 highest degree terms of the weighted Ehrhart 
polynomial for the weight h(x) = -^^^j— , we need only the fco+l lowest 
degree homogeneous terms of the meromorphic function S{Cs){^), for 
every vertex s of p. We compute such an approximation in the next 
section. 

4. Approximation of the generating function of a 

SIMPLICIAL AFFINE CONE 

Let c C be a simplicial cone with integral generators {vj,j = 
1, . . . ,d), and let s G Vq. Let ko < d. In this section we will obtain 
an expression for the /cq + 1 lowest degree homogeneous terms of the 
meromorphic function S{s + c)(^). Recall that if c is unimodular, the 
function S{s + c)(^) has a "short" expression, given by ([2]), 

S(.+ c)«) = e«'->n^-i^, 

i=i 

where {vi,i = 1, . . . ,d) are the primitive integral generators of the edges 
and s is the unique lattice point in the corresponding cell s + B. Thus 
in the unimodular case, computing the lowest degree components is 
immediate. 

When c is not unimodular, it is not possible to compute efficiently 
the Taylor expansion of the function S{s + B)(,^) at order ko, if the 
order is part of the input as well as the dimension d. In contrast, if 
the order k^ is fixed, we are going to obtain an expression for the 
terms of degree < ko which involves discrete summation over cones in 
dimension < k^ only, and determinants. For example, for ko = 0, the 
constant term S{s + B)(0) is the number of lattice points in the cell, 
which is equal to a determinant, by Lemma [6] (iii). 

We need some notations. 

For / C {1, . . . , d}, we denote by Lj the linear span of the vectors 
{vi,i G /). We denote by B/ = ^jg/[0, l[vi the unit cell in Lj. 
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We denote by I'^ the complement of / in {1, . . . , d}. We have V = 
Lj © Ljc. For X E V we denote the components by 

X = Xj + Xjc. 

Thus we identify the quotient V/Ljc with Lj and we denote the pro- 
jected lattice by A/ C L/. Note that L/ n A C A/, but the inclusion is 
strict in general. 

Example 12. vi = (1,0), fa = (1,2). Take I = {1}. Then A/ = Zf . 

We denote by C/ the projection of the cone c on the space Lj. Its 
edges are generated by Vj, j G /, and the corresponding unit cell B/ 
is the projection of B. Remark that Vj may be non primitive for the 
projected lattice A/, even if it is primitive for A, as we see in the 
previous example. This is the reason why in Lemma [6], we did not make 
the (unnecessary) assumption that the generators Vj are primitive. 

For u = (mi, . . . , Ud), we denote uj = ""j- 

We denote the binomial coefficient f^n^pij^y^ by (^) . 

Definition 13. Given a function I \—>- X{I) on the set of subsets I C 
{1, . . . ,d} with cardinal \I\ < ko, we denote 

T{s, c, fco, A)(0 = ^(^) ^o\{Bi.)S{sj + c,)(0(-l)'"l'l n 77^- 

Remark 14. The function S{si + Cj){C,) is a meromorphic function on 
the space L}. We extend it to V* by the decomposition V = Lj Q) Ljc 

It is easy to see that the function T{s, c, ko, X){^) lies in J^e{V*), 
and its expansion into homogeneous components has lowest degree —d. 
Thus 

T(s, c, ko, Xm = T[-d]is, c, ko, A)(0 + T[_,+i](s, c, ko, A)(0 + ■ ■ ■ • 

We will use functions I ^ X{I) which have the following property. 

Definition 15. A {d, ko)-patchfunction is a function I t— > A(/) on the 
set of subsets I C {l,...,d} of cardinal \I\ < ko which satisfies the 
following condition. 

(8) g«i+-+«d = ^ A(/)e"^ mod terms of u- degree >ko + l. 

\I\<ko 
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Theorem 16. Let 1 1— > A(/) be a ko-patchf unction. Then we have 
(9) Sis + 'Bm^Y. A(/)vol(B,.)5(s, + B,)(e) 

|/|<fco 

mod terms of S,- degree > /cq + 1- 

(10) 

S'(s + c)(,^) = T(s, c, fco) mod terms of ^-degree > —d + ko + l. 
Proof. Using ([2]), we write 

S{s + c)(0 =\S{s + B)(0 n 7^^^ ^ / , 

Thus we need only the terms of ^-degree at most in the Taylor expan- 
sion of the holomorphic function S{s + B)(^) 11^=1 finally 
we need only the the terms of ^-degree at most k^ in the Taylor expan- 
sion of S{s + B)(^). Applying ([8]) to each term e^^'^^ = e^^^'^^'"^^'^^'^ of 
the finite sum S{s + B)(^), we have 

5(s+B)(0=^ A(/) ^ e<«'^^> mod terms of ^degree > A;o+l. 

|/|<fco xe(s+B)nA 

For each /, the term X]2:e(s+B)nA e'^^'^^^ is the sum, over x G (s + B) flA, 
of a function of x which depends only on the projection xi of x in the 
decomposition x = xi + xjc & (B L^'^ . When x runs over (s + B) fl A, 
its projection xj runs over (s/ + B/) fl Aj. Let us show that the 
fibers have the same number of points, equal to vol(B/c). For a given 
xj G (s/ + B/) n A/, let us compute the fiber 

{ye (s + B) nA; yi = xi}. 

Fix a point xj + xic in this fiber. Then y = xj + yjc lies in the fiber 
if and only if yic — xjo G {sjc — xjo + B/c) fi A. By Lemma [6]^ii) , the 
cardinal of the fiber is equal to vol(B/c). Thus, we have obtained ([9]). 

Next we write the sum S{si + C/)(^) over the projected cone sj + Cj 
in terms of the sum over the projected cell sj + B/. We obtain 

S{s + cm^ E A(/)vol(B,Og(3/ + c,)(On n - e(g.-.)) 

^ J2 m^oi{Bj.)s{sj+cjm{-ir-^'^U7r—x 

\I\<ko jel^^^'^^' 

mod terms of ^-degree > —d + ko + 1. 
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□ 

Next we compute an explicit {d, /co)-patchfunction. 
Lemma 17. If I ^ {1, . . . ,d} has cardinal \I\ < k^, let 



1 1 

Let us consider P{t) := Hill + ^(^"' ~ 1)) = Sg=o ^'^(m)^^ as a polyno- 



mial in the indeterminate t. As e"' — 1 is a sum of terms of Wj-degree 
> 0, we have 

(11) gMi+- +Md = ^ Cq{u) mod terms of u-degree > fco + 1- 

g=0 

In order to compute the coefficient Cq{u), we write 




Then Xd,ko satisfies Condition 

Proof. The trick is to write e" = 1 + t{e^ — l)\t=i- Thus 



d 



d 



d d 



Pit) = Uil + t{e^' - 1)) = n((l 



1 1 



By expanding the product, we obtain 




Summing up these coefficients for < g < /cq, we obtain 




For mo < do, let us denote 




j=0 



Thus 



5^C,(n)= J2 Fiko-\Ild-\I\)e 



9=0 |/|<fco 
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The sum F{mo, do) is easy to compute by induction on mo, using the 
recursion relation 



We obtain 



F{mo,do) = (-i; 



mo 



do - 1 
mo 



Hence, 

The claim follows now from Equation ( |TT|) . □ 

Remark 18. As promised, the main feature of Formula fW\) is that the 
right-hand-side T(s, c, ko, A) involves discrete summations in dimension 
\I\ < ko only. 



Theorem [16] can be reformulated in terms of the intermediate valu- 
ations introduced by Barvinok in [3]. The reformulation relies on the 
next lemma, which shows that the [d, fco)-patchfunction condition is 
equivalent to a Moebius-type condition for the function / A(/). 

Lemma 19. Let < ko < d be two integers. Let X be a function on 
the set of subsets / C {1, . . . , rf} of cardinal \I\ < ko. The following 
conditions are equivalent, 
(i) For every Iq of cardinal |/o| ^ ^o? 

E ^(^) = 1- 

{/; |/|<fco, loCI} 

(a) For every integer k such that < k < ko, we have the equality of 
polynomials 

(12) (u, + --- + u,f= J2 ^W"/ 

|7|<fco 

(Hi) The function X is a {d, ko)-patchf unction. 

Proof. Conditions (ii) and (iii) are clearly equivalent. Let us show that 
(i) and (ii) are equivalent. We expand (mi + ■ ■ ■ + Ud)'' into a sum of 
monomials. A monomial of degree k can involve at most k variables 
Ui, with k < ko. Therefore we obtain 

(13) + E E Ili- 

\I\<ko {k,),ei iel * 
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We expand similarly each term in the right-hand side of f[T^ . A given 

monomial Hig/o 'kli with /Cj 7^ for all i ^ Iq, occurs in the right-hand 
side of f|T2l) exactly for the subsets I such that Iq C /. Thus (i) implies 
(ii). Conversely, Equation (|T2l) for k = ko implies (i). □ 

5. Computation of Ehrhart quasi-polynomials 

We now apply the approximation of the generating functions of the 
cones at vertices to the computation of the highest coefficients for a 
weighted Ehrhart polynomial, when the weight is a power of a linear 
form, as we explained in section [31 

Corollary 20. Let p be a simple lattice polytope. Fix ,^ G M"' and 
M eN. Let Em{p, ^,M),m = 0, . . . ,d + M, be the coefficients of the 
weighted Ehrhart polynomial 

xgnpnA m=0 

Fix < ko < d. Let X be a {d,ko)-patchf unction. Then, for k = 
0, . . . , /co; the Ehrhart coefficient EM+d-k{p,(,, M) is given by the fol- 
lowing formula. 

It: „\AI+d-k 

(14) i?M+d-.(p,e,M) = Jil^__^T(0,c„A:o, A) (0- 

In the general case, when the vertices are not assumed to be lattice 
points, we state the result without going through the details of the 
computation. 

Theorem 21. Let p be a simple polytope. For each vertex s of p, let 
Qs E N be the smallest integer such that qsS G A. For n G N, let rig be 
the residue of n mod qs, so that < Ug < Qs — 1. Fix ^ E V* and M a 
nonnegative integer. Fix < ko < d. Let X be a {d, ko)-patchf unction. 
Then the Ehrhart quasi-polynomial 

^ Ml 

x£npr\A 

coincides in degree > M + d — ko with the following quasi-polynomial 

It: „\M+d~k 

fc=0sev(p) ^ ^' 
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Observe that ( |T5l) is clearly a quasi-polynomial in n with coefficients 
which depend only on the residues (n mod g^) = ra^, s G V(p). 

Remark 22. In practice, we first reduce the vertices s mod A hy using 

S{s + c,)(0 = e^^^^^S{s -v + c,)(0, for v e A. 

Then we write an approximation similar to l[TS\) . 

6. Conclusion 

Let p C M'^ be a rational simplex. Let x) be a rational linear form 
on M.'^, and consider a power x)^. Let Em{p, M), m = 0, . . . , (i + 
M, be the coefficients of the weighted Ehrhart quasi-polynomial 

xgnpnA m=0 

Fix an integer k^, Q < < d. The main consequence of this study is a 
method for efficiently computing the + 1 highest degree coefficients 
Em{p^ ^5 M), for m = M + d, . . . , M + d — k^. The method relies on 
expanding ( fTSi) in Theorem [2T] as a power series in ^. 

Furthermore, one can write any homogeneous polynomial weight 
h{x) as a linear combination of powers of linear forms, 

h{x) = Y,cu{ik,xt'. 

k 

In a forthcoming article by the authors of [T], we will show how to 
derive 

- ffist, an algorithm for computing Em{p,C,,M), for m = M + 
d, . . . , M + d — ko- Hopefully this algorithm is polynomial, when the in- 
put consists of the dimension d and the degree M , the rational simplex 
p C M'^ , the rational linear form ^ on M*^, provided ko is fixed; 

- second, an algorithm for computing the feg + 1 highest degree coef- 
ficients of a weighted Ehrhart quasi-polynomial relative to a simplex. 
Hopefully this algorithm is polynomial when k^ and the degree of the 
weight h{x) are fixed, but the dimension of the simplex is part of the 
input. 

[1] dealt with the case of the highest Ehrhart coefficient which is just 
the integral of the weight over the simplex. 

References 

[1] Baldoni V., Berline N., De Loera J., Koppe M., Vergne M. How to 

Integrate a Polynomial over a Simplex. arXiv:0809.2083 



16 



VELLEDA BALDONI, NICOLE BERLINE, AND MICHELE VERGNE 



Baldoni V., Berline N. and Vergne M. Euler-Maclaurin expansion of 
Barvinok valuations and Ehrhart coefficients of rational polytopes 
Barvinok A. I., Computing the Ehrhart polynomial of a convex lattice poly- 
tope, Discrete Comput. geom. 12 (1994), 35-48. 

Barvinok A. 1., Computing the Ehrhart quasi-polynomial of a rational sim- 
plex, Mathematics of Computation, 75 (2006), 1449-1466. 
Barvinok A. I. and Pommersheim J., An algorithmic theory of lattice 
points in polyhedra, New Perspectives in Algebraic Combinatorics (Berke- 
ley,CA, 1996-97), Math. Sci. Res. Inst. Pubhc 38, Cambridge University Press, 
Cambridge, (1999), pp 91-147. 

Berline N. and Vergne M. Local Euler-Maclaurin formula for polytopes 
(2005), arXiv:math CO/0507256. To appear in Moscow Math. J. 
Brion M., Points entiers dans les polyedres convexes, Ann. Sci. Ecole Norm. 
Sup. 21 (1988), 653-663. 

De Loera J. A., Haws D., Hemmecke R., Huggins H., Tauzer J. and 
Yoshida R., A Users Cuide for LattE vl.l, 2003, software package LattE, 



available at http://www.math.ucdavis.edu~latte| 



Stanley R. Enumerative combinatorics. Vol 1 (1997), Cambridge Studies in 
Advanced Math. 49. 



Universita di Roma Tor Vergata, Dipartimento di Matematica, via 
DELLA Ricerca Scientifica, 00133 Roma, Italy 
E-mail address: baldoni@mat.uniroma2.it 

Ecole Polytechnique, Centre de mathematiques Laurent Schwartz, 
91128, Palaiseau, France 

E-mail address: berline@math.polytechnique.fr 

Institut de Mathematiques de Jussieu, Theorie des Groupes, Case 
7012, 2 Place Jussieu, 75251 Paris Cedex 05, France 

Ecole Polytechnique, Centre de mathematiques Laurent Schwartz, 
91128, Palaiseau, France 

E-mail address: vergne@math.polytechnique.fr 



